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ABSTRACT 

This study deals with the computational analysis of buoyancy-induced instability in the near- 
field of an isothennal helium jet injected into quiescent ambient air environment. Laminar, 
axisymmetric, unsteady flow conditions were considered for the analysis. The transport 
equations of helium mass fraction coupled with the conservation equations of mixture mass and 
momentum were solved using a staggered grid finite volume method. The jet Richardson 
numbers of 1.5 and 0.018 were considered to encompass both buoyant and inertial jet flow 
regimes. Buoyancy effects were isolated by initiating computations in Earth gravity and 
subsequently, reducing gravity to simulate the microgravity conditions. Computed results concur 
with experimental observations that the periodic flow oscillations observed in Earth gravity 
subside in microgravity. 
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1. INTRODUCTION 


Low-density gas jets are characterized by the injection of a lighter fluid into a dense ambient 
environment. The resulting flow structure is susceptible to buoyancy arising because of the 
difference in densities of jet and surrounding fluids. In general, jets are considered buoyant for 

Richardson number, Ri= gd(Poo — P-)/PjU. >1.0 and inertial for Ri«1.0, where g is the 

grav.tat.onal acceleration, d is the injector tube tnside diameter. p„ and p . are, respec.tvely, the 

ambient and jet fluid densities, and Uj is the average jet exit velocity. The present study 
examines low-density gas jets in both buoyant and inertial flow regimes, with the goal of 
delineating gravitational effects on the flow structure in the near-field. 

Several studies 1 ' 6 have focused on the far-field behavior of low-density jets, which is 
influenced by the instabilities arising in the near-field. Low-density gas jets are known to exhibit 
self-excited periodic oscillations in both buoyant and the inertial flow regimes. In the buoyant 
flow regime, Subbarao and Cantwell 7 conducted experiments with helium jets injected into a co- 
flow of air for Ri varying from 0.5 to 6.0. The frequency (f) of the periodic oscillations expressed 
by the Strouhal number, St=fd/Uj correlated with Richardson number for Ri>1.0 indicating 

o 

buoyancy dependent instability. Hamins et al. observed periodic oscillations in a helium jet 
injected into air, and reported that a minimum jet exit velocity was required to initiate the 
oscillations. Similar experiments were conducted by Cetegen and co-workers 9 ' 11 for 
axisymmetric and planar plumes, whereby buoyancy induced toroidal vortical structures 
contaminating the primary jet flow were observed. Diagnostics involved flow measurements 
using laser Doppler velocimetry and particle image velocimetry. Pasumarthi and Agrawal “ 
conducted experiments to characterize jet flow in terms of concentration measurements across 
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the whole field of a helium jet injected into quiescent air. Experiments by Yep et al. 13 
demonstrated that the self-excited flow oscillations in Earth gravity were absent in the 
microgravity environment of the 2.2s drop tower. 

Self-excited oscillations have also been reported in low-density inertial jets at Ri«1.0 14 " 17 . 
In these studies, buoyancy effects were considered negligible because of the small Richardson 
number. Although this assumption is justified in the jet core, buoyancy may be significant in the 
jet shear layer where density and velocity gradients are large. Recently, Pasumarthi 18 conducted 
experiments at Ri<1.0 using helium jet injected into quiescent air in the microgravity 
environment of the 2.2s drop tower. Concentration measurements revealed global oscillations in 
the near-field in Earth gravity. In microgravity, the oscillation amplitude decreased significantly 
for Ri<0.01 and no oscillations were observed for Ri>0.1. These experiments provided direct 
physical evidence for the first time that the flow oscillations in inertial low-density gas jets are 
buoyancy-induced, a well known phenomenon in buoyant jets. 

Past experiments have characterized buoyancy effects using either concentration or velocity 
measurements. However, simultaneous visualization of velocity and concentration fields is 
desired to fully understand the dynamics of the oscillating flow structure. Experimental 
difficulties in obtaining simultaneous vector-scalar measurements, especially in the low-gravity 
environment, have prompted studies of gas jets using computational fluid dynamic (CFD) 
analysis. Mell et al. 19 performed numerical simulations of a helium jet injected into air and 
found that the computed oscillation frequency matched with experiments of Hamins et al . 
Recently, Soteriou et al2 investigated the near-field dynamics of planar buoyant plumes 
incorporating a Lagrangian transport element method. Computations qualitatively captured the 
plume instantaneous behavior observed experimentally. The mechanism of vorticity generation 
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and relative roles played by buoyancy and viscous forces in explaining the plume instability were 
discussed. However, a description of the flow field and its interaction with the scalar structure is 
desired to explain the instability in low-density gas jets. On a related note, several numerical 
studies have been conducted to understand the dynamics of flickering jet diffusion flames" ’ . 
Thus, the objective of the present study is to delineate the effects of buoyancy in low-density gas 
jets. The CFD analysis is utilized to concurrently visualize the velocity and concentration fields. 
Both buoyant and inertial jets are examined in Earth gravity, microgravity, and during change 
from Earth to microgravity. 


2, NUMERICAL FORMULATION 


2. 1 Governing Equations 

Consistent with experimental studies ’ , the numerical formulation is based on unsteady, 

laminar, and axisymmetric flow behavior. The following additional approximations were made 
to simplify the analysis: (i) the flow is incompressible, (ii) temperature and pressure variations 
are negligible and (iii) the flow is a binary fluid system with air treated as a single species. 
Accordingly, the governing equations are expressed as follows: 


dp 15/ x d / x 

¥ + 7 ar(^) + &(H=o 
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Equation ( 1 ) represents the continuity equation and Equation (2) is the generalized transport 
equation, which represents the momentum, or species conservation equation depending on the 
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variable F 1 ’ and S* are, respectively, the transport coefficient and source term of the 
variable^ , as listed in Table 1. 


Table 1. Transport coefficients and source terms in Equation 2 
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Density was computed from the equation of state for incompressible ideal gas. The dynamic 

23 

viscosity was calculated using Wilke’s mixture averaged formula” . The binary diffusion 
coefficient of helium into air, Db was specified as a constant. In microgravity, the gravitational 
acceleration was specified as lOOpg, where g is the gravitational acceleration on Earth. 

2.2 Boundary Conditions 

The computational domain extended 15d in the axial direction and 6d in the radial direction 
where “d” is the tube inside diameter. These dimensions were established after trials indicating 
that the near-field flow behavior was independent of the domain size. The upstream boundary 
was placed at a distance of Id from the jet exit plane to account for the diffusion upstream of the 
jet exit observed in experimental studies . Symmetric boundary condition was imposed along 
the jet centerline. At the tube inlet, the velocity profile was fully developed and the helium mole 
fraction was unity. No slip and zero mass diffusion conditions were imposed at the tube walls. 
At the exit boundary, pressure outlet condition was implemented where flow properties are 


5 




extrapolated from the interior 24 . Pressure inlet and outlet conditions were specified, respectively, 
at the upstream and far-field radial boundaries. All flow properties were set to zero at the start of 
the computations. The operating pressure and temperature were specified, respectively, as 1 atm 
and 300K. 

2.3 Computational Procedure 

An orthogonal non-uniform grid system was used to split the computational domain into five 
sub regions. The grid elements in the radial direction were concentrated in the flow oscillation 
region. The grid size was increased gradually in the axial direction. The governing equations 
were solved sequentially using a segregated approach 24 . The pressure -velocity coupling was 
achieved using the SIMPLEC algorithm. Grid sensitivity analysis was performed to obtain grid- 
independent solution. Striking a balance between computational time and accuracy, a grid with 
31093 and 80021 nodes was used, respectively, for buoyant and inertial jet analysis. The time 
step for buoyant and inertial jet computations was 1.6ms and 0.3ms, respectively. The small time 
step was needed to resolve the high-frequency oscillations in the inertial jet. The numerical code 
was validated with experimental data ’ and experimental-numerical comparisons are presented 
for Earth gravity conditions in Ref. 26 and for change from Earth gravity to microgravity in later 
sections of this paper. 


3. RESULTS AND DISCUSSION 

Time-dependent simulations of the isothennal helium jet injected into the ambient air 
environment were performed with the objective of quantifying buoyancy effects in the near-field. 
Although computations were performed for several cases, detailed results are presented for only 
one buoyant jet and one inertial jet at operating conditions listed in Table 1. First, details 
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pertaining to the flow behavior of a buoyant helium jet in Earth gravity and during change from 
Earth gravity to microgravity are presented. Similar information is provided for the inertial jet. 
Finally, flow behavior of buoyant and inertial jets is compared in Earth gravity and microgravity. 


Table 1. Summary of Test Conditions 


Case 

Re 

Ri 

d, mm 

Uj, m/s 

f > Hz 
(Computed) 

fHz 

(Measured) 

Case 1 (Buoyant) 

300 

1.52 

31.8 

1.15 

14.5 

12.5 

Case 2(Inertial) 

800 

0.018 

14.05 

6.75 

150 

120 


3.1 Buoyant Jet 

3.1.1 Flow and Concentration Fields in Earth Gravity 

Figure 1 shows a sequence of helium mole percentage contours superimposed with 
velocity vectors for the buoyant jet (case 1). At t=0.0ms, buoyancy accelerates the jet core, 
which contracts to conserve the mass. This results in entrainment of ambient fluid producing a 
toroidal vortex represented by the flow recirculation. Note that the vortex core represented by a 
black dot is located at z/d=0.12 in Fig.l. In the next plot at t=13.6ms, the vortex is characterized 
by a larger recirculation region accompanied with greater contraction of the jet core, evident 
from the inward indentations in the concentration field near z/d=0.2. Subsequently, between 
t=27.3ms and 54.7ms, the vortex grows in size and convects downstream to contaminate a 
greater portion of the jet core. The vortex gains its strength by momentum transfer from the jet 
core. This feature is evident in Fig. 1 wherein the velocity vectors downstream of the vortex core 
are larger than those upstream before interacting with the jet core. As the vortex propagates 
downstream, the jet expands near the tube exit inducting buoyant fluid to initiate another vortex 
at t=68.4ms to repeat the oscillation cycle. These results illustrate self-excited periodicity in the 
jet at a frequency of 14.6Hz, in excellent agreement with the oscillation frequency of 12.5Hz 

13 

measured by Yep et al. 
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3.1.2 Flow and Concentration Fields during Change from Earth gravity to Microgravity 

Computations were perfonned to visualize the jet flow during change from Earth gravity to 
microgravity and to depict how the jet flow adjusts itself in the absence of buoyancy. As a first 
step, the simulation results were compared with experiments to validate the computational 
model. Figure 2 shows a series of computational-experimental helium mole fraction contour 
plots during change from Earth gravity to microgravity. At the onset of microgravity denoted as 
T=0.0s, the indentations in the helium concentration level higher than 50% signify a vortex 
located at z/d=0.7. After change to microgravity at T=16.66ms, the vortex convects downstream 
as the jet expands in the near exit region. At T=50.0ms, the jet has expanded at the exit to r/d= 
0.7. The indentations in the concentration contours have weakened, signifying the diminishing 
nature of the vortex. Subsequently, the jet widens gradually throughout the near-field with 
diminishing oscillations as steady conditions are reached in microgravity. The helium mole 
fraction contours in microgravity are straight lines reminiscent of a non-buoyant jet. The above- 
mentioned features replicate similar phenomena observed in experiments showing self-excited 
oscillations gradually subsiding after microgravity was initiated. Judging from the plots in Figure 
2, it is inferred that computations reproduced experimental observations, proving model’s 
capability in predicting the flow structure during change from Earth gravity to microgravity. 

Figure 3 shows contours of helium mole percentage overlapped by velocity vectors during 
change from Earth gravity to microgravity to highlight the interaction of the flow field with the 
concentration field. At T=0.0ms, a vortex with core located at z/d=0.7 characterizes the 
entrainment of ambient fluid into the jet. At this instant the jet flow is similar to that at 
T=27.3ms in Fig. 1(c). The jet width characterized by the 10% helium concentration contour 
extends to r/d=0.6 near the exit. In microgravity at T=33.33ms, the jet at the exit has expanded to 
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r/d=0.7. The vortex has convected downstream with its core located at z/d=1.8. Velocity vectors 
show a weakening of the vortical structure at this instant. At T=66.6ms, the weakened vortex has 
moved downstream of the field-of-view while another vortex is fonned near the jet exit at 
z/d=0.3. At T=133.3ms, flow recirculation has diminished as the jet expands to r/d=0.78 near the 
exit. Similar features are exhibited at t= 166.6ms with the flow approaching steady conditions in 
microgravity. Finally, at T=2000ms the jet has assumed a steady columnar shape. The flow 
velocity decreases in the radial direction as the jet fluid mixes with ambient air. Further details 
about the flow structure are provided by spatio-temporal plots during the change from Earth 
gravity to microgravity. 

3.1.3 Temporal Evolution of Axial Velocity and Concentration Fields 

Figure 4 depicts temporal evolution of the concentration field during the change from Earth 
gravity to microgravity. Computed and measured 1 3 time traces of helium mole fraction at z/d=1.0 
illustrate that the oscillating jet flow expands radially after the change in gravity. Results show 
that the oscillations sustain during a brief initial period of microgravity, with steady conditions 
reached within T=0.75s. Again, the experimental trends are predicted well by the computations. 
Figure 5 shows time trace plots of helium mole percentage and axial velocity at selected 
locations during the change from Earth gravity to microgravity. Near the jet exit at z/d=0.15, the 
helium concentration at the jet center (r/d=0) is constant at 100% throughout. In contrast, the 
axial velocity at the jet center fluctuates between 2.3m/s and 2.6m/s in Earth gravity. In 
microgravity, the axial velocity initially decreases and then reaches a steady value of 2.2m/s at 
T=0.1s. In the entrainment region at r/d=0.45, the helium mole percentage in Earth gravity 
varied between 45% and 94%. During change to microgravity, the helium concentration peaked 
to 100% within 0.08s. Afterwards, minor oscillations occurred before a steady value of 95% was 
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reached within T=0.3s. Correspondingly, the axial velocity reached a steady value of 0.4m/s in 
microgravity. Near the jet boundary at r/d=0.6, low amplitude oscillations in the concentration 
field are observed in Earth gravity. In microgravity, the helium mole percentage rose to 70% 
prior to reaching a steady value of 52% within T=1.0s. In Earth gravity, the axial velocity profile 
at this radial location shows troughs with negative velocity indicating flow recirculation. 

At a downstream location, z/d=2.0, pure helium is present throughout at the jet center. The 
axial velocity and its oscillation amplitude have increased at this location, signifying buoyant 
acceleration in Earth gravity. The flow oscillations diminish within 0.5s of microgravity as 
steady value of 2.1m/s is reached. At r/d=0.45, high amplitude concentration oscillations are 
observed in Earth gravity. At r/d=0.6, helium reaches peak concentration level of 60% during 
change from Earth gravity to microgravity and reaches a steady value of 56% at T=0.6s. The 
axial velocity profiles at the above radial locations depict flow behavior similar to that observed 
upstream. These features agree with time traces of concentration field measured by Yep et al. 13 
3.2 Inertial Jet 

3.2.1 Flow and Concentration Fields in Earth Gravity 

Figure 6 shows a sequence of helium mole percentage contours superimposed with 
velocity vectors to depict an oscillation cycle. In this case, the computed results are presented 
only between r/d=0.3 and 0.65 to visualize the flow oscillation region in greater detail. At t=0.0s, 
buoyant acceleration contracts the jet core leading to entrainment and formation of a toroidal 
vortex. The flow structure during the oscillation cycle is similar to that of the buoyant jet. The 
nature and propagation of the toroidal vortex are the same as already explained. 

Although self-excited flow oscillations were observed in both buoyant and inertial jets, some 
aspects of the flow field were found to be different. First, based on the concentration field, the 
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inertial jet is wider compared to the buoyant jet. This implies that helium diffuses farther in an 
inertial jet. This effect is opposite to constant density jets where the jet width decreases with 
increasing initial jet momentum. Contraction and mixing by entrainment are confined to a 
narrow region in the inertial case compared to the buoyant case. The oscillation frequency for the 
inertial jet (150 Hz) is much higher than that for the buoyant case (14.5Hz). The above results 
illustrate self-excited periodicity in the flow field irrespective of the initial jet momentum, which 
agrees qualitatively and quantitatively with experiments 13 ' 18 . 

3.2.2 Flow and Concentration Fields during change from Earth gravity to Microgravity 
Figure 7 shows the contours of helium mole percentage overlapped by velocity vectors 

during change from Earth gravity to microgravity and then to zero gravity to highlight the 

interaction of flow field with concentration field. At the onset of microgravity (T=0.0ms), a 

toroidal vortex is located at z/d= 0.75. The vortex convects downstream and a new one is formed 

near the exit as shown in Fig. 7(b) for T=24ms. Subsequently, the jet flow maintains a periodic 

structure in microgravity as by Figs. 7(c)-(e) representing an oscillation cycle. The curvature of 

concentration contours is less severe and the vortex is narrower in microgravity compared to 

.18 

those in Earth gravity. This result is consistent with experiments by Pasumarthi showing flow 
oscillations subsiding to smaller amplitudes in microgravity for Ri<0.01. The last plot in Fig 7 
pertains to the results obtained in zero gravity, simulated by completely turning off the gravity. 
In this case, the jet assumes a steady non-oscillatory state to support buoyancy as the cause of the 
instability. Results show that minor buoyant acceleration destabilizes the inertial jet. 

3.2.3 Temporal Evolution of Axial Velocity and Concentration Fields 

Figure 8 shows time trace plots of helium mole percentage and axial velocity at selected 
locations during change from Earth gravity to microgravity. At z/d=0.5, helium concentration at 
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the jet centerline remained 100% throughout. Concentration fluctuations for the inertial case are 
smaller compared to the buoyant case. Axial velocity at the jet center oscillates between 13 and 
14 m/s in Earth gravity and between 12.5 and 13.5 m/s in microgravity. These oscillations ceased 
when the gravity was completely turned off in zero gravity. At the downstream location, z/d=2.0, 
the axial velocity fluctuates between 12 and 15 m/s in Earth gravity. These fluctuations are 
slightly smaller in microgravity but absent in zero gravity. The concentration and velocity fields 
in Fig. 8 reveal that the flow oscillations subside in microgravity and disappear in zero gravity. 
These results ascertain that buoyancy plays an important role in both buoyant and inertial low- 
density gas jets. 

3.3 Comparison of Flow Behavior in Buoyant and Inertial Jets 

In this section, buoyant and inertial jet flows are compared to identify similarities and 
differences in the flow behavior. 

3.3.1 Vortex Convection Velocity 

Figure 9 shows the vortex convection velocity (U v ) nonnalized by the averaged jet exit 
velocity in Earth gravity. This plot was generated by tracking the position of the vortex core at 
different times. Results show a gradual increase of vortex convection velocity in the flow 
direction. The plots are similar for buoyant and inertial jets although higher values are reached 
for the buoyant case. 

3.3.2 Axial Velocity Profiles 

As already stated, buoyancy plays an important role in describing the instability phenomena. 
In case of the buoyant jet, buoyancy effects are prominent in the jet core. However, in the inertial 
jet, buoyancy affects the low-momentum region of the jet shear layer. This phenomenon is 
illustrated by the profiles of instantaneous axial velocity at various radial locations as shown in 


12 



Fig. 10. The axial velocity was normalized by the jet exit velocity (U e ) at that particular radial 
location. The profiles show the axial velocity varying with axial coordinate and time 
corresponding to the phases in Figs. 1 and 6. 

The normalized axial velocity profiles for the buoyant case (Figs. lOa-c) show periods of 
acceleration and deceleration at any given phase. The phases reveal the periodic flow behavior of 
the near-field region. Along the jet centerline, at phase A, slight deceleration occurs up to 
z/d=0.8 followed by a rapid acceleration between z/d=1.0 and 2.0. Note that the velocity peak at 
z/d=2.0 occurs slightly upstream of the vortex core at z/d=2.15. Near the jet exit, the decelerating 
flow at phase A accelerates at Phase B as a new vortex is formed. At phase C, rapid acceleration 
is observed near the exit, with peak normalized axial velocity reaching a value of 1.25 at z/d=0.5 
or slightly upstream of the vortex core at z/d=0.80. Subsequently, the flow decelerates because of 
the transfer of the energy from the jet core to the vortex. At later phases (D and E), the flow 
decelerates near the exit but accelerates at downstream locations. Note that the normalized peak 
axial velocity has increased to about 1.7 at phase D and to 2.0 at phase E. A common feature is 
that the flow accelerates upstream of the vortex core and decelerates downstream of it. The 
inertial jet shows a similar behavior. However, the normalized axial velocity profiles are 
elongated implying that the peak occurs farther downstream. The peak nonnalized axial velocity 
is within 1.1 for all phases. These results clearly show that the buoyancy induced acceleration at 
the jet center is more prominent for the buoyant case compared to the inertial case. 

At r/d=0.3, both buoyant and inertial jets exhibit flow acceleration and deceleration because 
of the interaction with the outer vortex. However, the peak nonnalized axial velocity has 
decreased from 2.0 to 1.4 for the buoyant jet and it has increased from 1.1 to 1.2 for the inertial 
jet. Thus buoyancy is more important at this radial location for the inertial jet. At r/d=0.45, both 
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buoyant and inertial jets show large variations in the normalized axial velocity because of the 
small initial momentum in the shear layer region. Note that the flow acceleration is greater for 
the inertial jet compared to that for the buoyant case. 

3.3.3 Mean and RMS Flow Structure 

The oscillating flow of buoyant and inertial jets is compared in Fig. 11 using the mean and 
RMS profiles of helium mole percentage and axial velocity at several axial locations. Near the jet 
exit at z/d=0.05, the mean concentration shows a top hat profile with finite shear layer thickness 
for both buoyant and inertial jets. Pure helium is present till r/d=0.38 for the buoyant case 
compared to r/d=0.46 for the inertial case. A steep reduction in helium concentration occurs 
thereafter. The scalar jet width extends to r/d=0.8 for the buoyant case and to r/d>0.8 for the 
inertial case, indicating greater radial expansion at the higher Reynolds number. 

At z/d=0.5, the mean concentration at the tube inside the boundary plane (r/d=0.5) has 
reduced to 50% for the buoyant case compared to 75% for the inertial case, suggesting greater 
radial mixing for the buoyant case. Further downstream at z/d=1.5, the concentration profiles 
show an inflection point located at r/d=0.29 for the buoyant case and r/d=0.33 for the inertial 
case, attributed to the entrainment by the vortex. The concentration gradient in the inflection 
region is higher for the buoyant case compared to the inertial case. The scalar jet width is nearly 
independent of the axial location for the buoyant case because the entrainment has restricted 
radial expansion of the jet. In contrast, for the inertial case, the scalar jet width increases in the 
flow direction because of the higher jet momentum limiting the entrainment. 

The RMS concentration profiles show large fluctuations near the jet exit for the buoyant 
case. The peak RMS concentration of 14% at z/d=0.05 increased to 25% at z/d=0.5 and to 29% 
at z/d=1.5. The location of the peak RMS concentration shifts inwards with the axial direction, 
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which is consistent with the contraction of the jet as the vortex propagates downstream. The 
concentration fluctuations confined to a wake region (0.4<r/d<0.6) near the jet exit broaden in 
the flow direction and contaminate a greater portion of the flow at z/d>1.0. These observations 
agree with experimental results presented by Pasumarthi. For the inertial case, fluctuations are 
negligible at z/d=0.05. Peak RMS concentration of 10% at z/d=0.5 increased to 12% at z/d=1.5. 
Lower fluctuations in the inertial jet compared with buoyant jet signify low level of mixing in the 
inertial jet. 

Mean axial velocity at z/d=0.05 shows a typical parabolic profile at the tube exit for buoyant 
and inertial cases. For the buoyant case, the axial velocity at the jet centerline has increased from 
2.4m/s at z/d=0.05 to 2.6 m/s at z/d=0.5 and to 3.1 m/s at z/d=1.5, indicating strong buoyancy 
induced acceleration. However, the buoyant acceleration is confined mainly to r/d<0.3. For the 
inertial jet case, the mean axial velocity is nearly the same at all axial locations. The RMS axial 
velocity profiles show negligible fluctuations near the jet exit (z/d=0.05) for both cases. The peak 
fluctuations increase with axial distance. For the buoyant case, peak RMS velocity of 0.3m/s at 
z/d=0.5 increased to 0.75m/s at z/d=1.5. Note that the peak RMS axial velocity shows a 
decreasing trend from the jet centerline towards the jet boundary plane for the buoyant case. In 
contrast, the RMS axial velocity peaked at r/d= 0.4 for the inertial case. 

3.3.4 Flow Behavior in Microgravity 

Radial profiles of axial velocity and concentration field in microgravity are presented in Fig. 
12 for buoyant and inertial cases. These profiles were obtained after the jet adapted to the change 
in gravity. In microgravity, the jet assumed steady, non-oscillatory structure for the buoyant case 
and it maintained flow oscillations for the inertial jet. For both cases, the jet width increased in 
microgravity. For buoyant jet, the helium concentration at r/d=0.7 is less than 5% in Earth 
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gravity and approximately 30% in microgravity. The corresponding values for the inertial jet are 
between 10-30% in Earth gravity and between 25-50% in microgravity. Figure 12 shows that the 
oscillation amplitude for the inertial jet is significantly smaller in microgravity compared to that 
in Earth gravity. 

CONCLUSIONS 

A laminar, time-dependent, axisymmetric numerical model was developed to simulate 
helium-air jets previously studied experimentally. The CFD-based model accurately simulated 
the experiments in Earth gravity and microgravity involving buoyant and inertial jets. The major 
conclusions of the study are summarized in the following: 

• The buoyant jet exhibited self-excited oscillations in Earth gravity at a frequency of 
14.5Hz. The inertial jet exhibited similar self-excited oscillations in Earth gravity at a 
frequency of 150Hz. However, the oscillation amplitude was smaller and oscillations 
were confined to a narrow region in the inertial jet compared to the buoyant jet. 

• The buoyant jet attained a steady, columnar structure in microgravity. Temporal plots of 
velocity and concentration fields depict that the flow oscillations persisted for a brief 
period during the change from Earth gravity to microgravity. 

• The inertial jet exhibited self-excited oscillations in microgravity, although the oscillation 
amplitude was lower. Simulations illustrated that the inertial jet reached a steady, non- 
oscillatory flow when the gravity was completely turned off (zero gravity). These results 
show that the inertial jet becomes unstable in the presence of small amount of buoyant 
acceleration. 
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• In Earth gravity, the scalar jet width was found to increase with increasing Reynolds 
number. This near-field behavior of low-density jets is opposite to that of constant 
density jets where the jet width decreases at higher Reynolds numbers. 

• The scalar jet width increased in microgravity for both buoyant and inertial jets. 
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Figure 1 . Contours of helium mole percentage superimposed by velocity vectors during an oscillation cycle in 
Earth gravity for Re=300, Ri=1.52, and d=31.8mm. (a) t=0.0ms; (b) t=13.6ms; (c) t=27.3ms; (d) t=41.0ms; (e) 
t=54.7ms; (f) t=68.4ms. 

Figure 2. Measured (1) and Computed (2) helium mole fraction contours during change from Earth gravity to 
microgravity. Re = 300, Ri = 1.52, and d=3 1.8mm. (a) T=0.0ms; (b) T=16.66ms; (c) T=50.00ms; (d) 
T=2 16.66ms; (e) T=2000ms. 

Figure 3. Contours of helium mole percentage superimposed by velocity vectors during change from Earth 
gravity to microgravity for Re=300, Ri=1.52, and d=3 1.8mm (a) T=0.0ms; (b) T=33.33ms; (c) T=66.66ms; (d) 
T=133.33ms; (e) T=166.66ms; (f) T=2000ms. 

Figure 4. Measured (a) and Computed (b) time traces of helium mole fraction during change from Earth gravity 
to microgravity. Re=300, Ri=1.52 and d=3 1.8mm at z/d=1.0. 

Figure 5. Time traces of 1) helium mole percentage and 2) axial velocity during change from Earth gravity to 
microgravity for Re=300, Ri=1.52 and d=31.8mm. a) z/d=0.15 b) z/d=2.0. 

Figure 6. Contours of helium mole percentage superimposed by velocity vectors during an oscillation cycle in 
Earth gravity for Re=800, Ri=0.018 and d=14.05mm. (a) t=0.0ms; (b) t=1.2ms; (c) t=2.4ms; (d) t=3.6ms; (e) 
t=4.8ms; (f) t=6.6ms. 

Figure 7. Contours of helium mole percentage superimposed by velocity vectors during change from Earth 
gravity to microgravity for Re=800, Ri=0.018 and d=14.05mm. (a) T=0.0ms; (b) T=24ms; (c) T=400ms; (d) 
T=402.5ms; (e) T=405.8ms; (f) zero gravity. 

Figure 8. Time traces of 1) helium mole percentage and 2) axial velocity during change from Earth gravity to 
micro gravity for Re=800, Ri=0.018 and d=14.05mm at a) z/d=0.5 b) z/d=2.0. 

Figure 9. Nonnalized Vortex Convection velocities for Re=300, Ri=1.52 and d=31.8mm and Re=800, Ri=0.018 
and d=l 4.05mm cases. 

Figure 10. Nonnalized axial velocity for 1) Re=300 and 2) Re=800 at different radial locations of a) r/d=0.0 b) 
r/d=0.3 c) r/d=0.45. 

Figure 11. Mean and RMS helium mole percentage and axial velocity profiles at different axial stations for a) 
Re=300, Ri=1.52, and d=3 1.8mm and b) Re=800, Ri=0.018, d= 14.05mm in Earth gravity. 

Figure 12. Helium mole percentage and axial velocity profiles at different axial stations for a) Re=300, Ri=1.52, 
d=3 1.8mm and b) Re=800, Ri=0.018, d=l 4.05mm cases in microgravity. 
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